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ABSTRACT 

Context. We present the results from a study of the X-ray variability and the near-IR to X-ray spectral energy 
distribution of four low-luminosity, Seyfert 1 galaxies. 

Aims. We compared their variability amplitude and broad band spectrum with those of more luminous AGN in order 
to investigate whether accretion in low-luminosity AGN operates as in their luminous counterparts. 
Methods. We used archival XMM-Newton and, in two cases, ASCA data to estimate their X-ray variability amplitude 
and determine their X-ray spectral shape and luminosity. We also used archival HST data to measure their optical 
nuclear luminosity, and near-IR measurements from the literature, in order to construct their near-IR to X-ray spectra. 
Results. The X-ray variability amplitude of the four Seyferts is what one would expect, given their black hole masses. 
Their near-IR to X-ray spectrum has the same shape as the spectrum of quasars which are 10 2 — 10 5 times more 
luminous. 

Conclusions. The objects in our sample are optically classified as Seyfert 1-1.5. This implies that they host a relatively 
unobscured AGN-like nucleus. They are also of low luminosity and accrete at a low rate. They are therefore good 
candidates to detect radiation from an inefficient accretion process. However, our results suggest that they are similar 
to AGN which are 10 2 - 10 s times more luminous. The combination of a "radiative efficient accretion disc plus an 
X-ray producing hot corona" may persist at low accretion rates as well. 
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qq ' 1. Introduction tion in these objects operates in a radiative inefficient mode. 

, For example, according to Merloni et. al. (2003), the accre- 

L| ■ The current paradigm for active galactic nuclei (AGN) com- tion mode should change below L x /L Ed d ~ 10~ 3 (where 

. £^ \ prises a supermassive black hole (BH) which is accreting L x and L E dd are the 2-10 keV and Eddington luminos- 

^ . material through a disc. Accretion in a geometrically thin, ity, respectively), while Chiaberge, Capetti & Macchetto 

' optically thick disc (Shakura & Sunyaev 1973) is the fa- (2005) suggest that RIAFs may operate in active nuclei 

Ctf \ vorite mechanism to achieve the high accretion efficiency (of with L /L E dd < 10 -4 (where L is the optical luminosity 

the order of ~ 0.1 or so) which is necessary to explain the in the R— band). However, conclusive evidence for the pres- 

large power emitted by AGN. Theoretical modeling of low ence of RIAFs in low luminosity AGNs (LLAGNs) is still 

accretion rate systems suggests that, when the accretion lacking. 

rate is below a few percent of the Eddington limit, L E dd, One of the most frequently applied tests to investigate 
AGN may switch to a different accretion mode character- whether a RIAF operates in LLAGNs is to compare their 
izcd by low radiative efficiency. Most "radiatively inefficient nuclear spectral energy distribution (SED) in the near- 
accretion flow" models (RIAFs; see, e.g., Yuan (2007) and infrared (IR) to ultraviolet (UV) part of the spectrum with 
references therein) suggest that the kinetic energy associ- the average SED of powerful quasars. At these frequencies 
ated with the gas is either advected with the matter into different accretion disc models are expected to show large 
the BH or redirected into an outflow. differences in spectral shape. RIAFs should lack both the 
Many of the nearby AGN are intrinsically faint (ex- "big blue bump" and the IR (reprocessed) bump, which in- 
hibiting bolometric luminosities less than io 43-44 ergs/s stead characterize the emission from an optically thick, ge- 
or so) while the central BH in a large fraction of them is ometrically thin accretion disc and the surrounding heated 
as large as 10 7 — 10 8 M© (see, e.g., Panessa et al. 2006). dust. Application of this test to various samples of LLAGNs 
Consequently, most of these nuclei accrete at a rate lower has yielded contradictory results. While Maoz (2007) finds 
than a few percent of L E dd- It is then possible that accre- that the broad band SEDs of AGNs with luminosities as low 

as ~ 10 40 erg/s are quite similar to the SED of more lumi- 

Send offprint requests to: I.E. Papadakis; e-mail: nous AGN, Ho (1999) found that the low-luminosity AGN 

jhep@physics.uoc.gr SEDs have a weak or absent blue bump and are "radio 
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loud". Quataert, di Matteo and Narayan (1999) find that 
the optical/UV spectrum of M81 and NGC 4579 decrease 
with increasing frequency, in contrast to the "canonical" 
quasar spectrum. Ptak et al. (2004) find that NGC 3998 (a 
"type 1" LINER galaxy) is lacking the "big blue bump", 
while RIAF models can fit its UV to X-ray SED reasonably 
well. Chiaberge et al. (2006) also find that the SED of a low 
luminosity Seyfert 2 galaxy, namely NGC 4565, is different 
from that of luminous quasars. 

Another powerful test to identify the presence of RIAFs 
in low luminosity AGNs is to compare their variability prop- 
erties with the properties of more luminous AGNs. Ptak et 
al. (1998) found that LLAGNs tend to show little or no 
significant short-term variability, despite their small lumi- 
nosity. This is opposite to what is observed in powerful 
AGN, where the variability amplitude increases with de- 
creasing luminosity. This difference has been interpreted as 
an indication that the size of the X-ray producing region 
in LLAGNs is significantly larger than the size of the X- 
ray source in "normal" Seyferts. This is in support of the 
hypothesis that advection-dominated accretion operates in 
LLAGNs, as in this case one can explain the larger size of 
the nuclear source. 

Our aim is to compare the X-ray variability amplitude 
and the optical/X-rays SED of low luminosity AGN and 
powerful quasars. We used the Palomar optical spectro- 
scopic survey of nearby galaxies (Ho, Filippenko & Sargent 
1995) to find suitable objects for our study. This sample has 
the advantage of having uniform and high-quality data that 
allowed classification of the galactic nuclei to be determined 
with well defined and objective criteria (Ho, Filippenko & 
Sargent 1997). We only considered objects which are clas- 
sified as Seyfert 1 - Seyfert 1.5, because in this case we can 
be certain of the AGN- like nature of their nuclear activity, 
and that the central source is unobscured. There are 11 such 
objects in the Palomar sample (Ho et al. 1997). We chose 
four from these sources, namely NGC 4235, NGC 4639, 
NGC 5033 and NGC 5273, because there exist estimates of 
their central black hole mass (hence we can estimate their 
LEdd) fmd their nuclear luminosity is small: their Lx/LEdd 
ratio is smaller than 10 -3 . Therefore, it is possible that ac- 
cretion in these objects operates in a radiatively inefficient 
mode. 

First, we used archival XMM-Newton and ASCA data to 
measure their X-ray variability amplitude and we compared 
it with the variability amplitude of more luminous Seyferts. 
Then, by i) fitting models to the XMM-Newton spectra of 
the four sources, ii) measuring their nuclear optical lumi- 
nosity using archival HST data, and iii) using near-IR lu- 
minosity measurements from the literature, we constructed 
their near-IR to X-ray SEDs and we compared them with 
the mean, broad band SED of luminous PG quasars. Our 
final result that the properties of the four objects in our 
sample are similar to those of more luminous AGNs, should 
have interesting implications regarding the nature of the ac- 
cretion mode in low luminosity active galaxies. 

2. The X-ray timing and spectral analysis 

2.1. The XMM-Newton observations 

All four Seyferts have been observed by XMM-Newton. 
Although the exposures are less than ~ 15 ksec long, XMM- 
Newton can provide high quality data for the detection of 



short term variations and the accurate determination of 
their nuclear, X-ray spectral shape and luminosity. 

We have retrieved data from the public XMM-Newton 
archive. Details of the observations are given in Table [TJ 
The EPIC data were reprocessed with the XMMSAS version 
6.5. None of the observations were affected by strong back- 
ground flaring events. With an average count rate of a few 
cts s _1 photon pile-up is negligible for the PN and the MOS 
detectors, as verified using the XMMSAS task epatplot. 

For the PN instrument, source counts were accumulated 
from a rectangular box of 27x26 RAW pixels (1 RAW pixel 
~ 4.1") around the position of the source. Background 
data were extracted from a similar, source free, region 
on the chip. We selected only single and double events 
(PATTERN<4, FLAG<0; for details of the instruments see 
Ehle et al. 2007) in the energy range from 400 eV to 10 
keV. For the MOS data we accumulated the source counts 
using a circular aperture with radius 60" centered on the 
source positions. The background data were extracted from 
a similar region on the same chip. Events with PATTERN 
< 12 and FLAG=0 were used for the analysis. 

The XMM-Newton images for all objects show a bright, 
compact, unresolved nucleus. NGC 4639, NGC 5033 and 
NGC 5273 have also been observed with Chandra. The 
higher resolution Chandra images reveal the presence of a 
bright nucleus in all cases as well. Weak diffuse emission, 
which is less prominent in the hard band, appears in the im- 
ages of NGC 4639 and NGC 5033. There are two off-nuclear 
sources (at a distance of less than ~ 1 arcmin) in the case 
of NGC 5033 and one such source in NGC 5273. These 
sources are significantly weaker than the central nucleus 
in both cases. We therefore expect that the XMM-Newton 
data should not be significantly affected by contamination 
from either unresolved sources or diffuse emission. 

2.2. Timing analysis 

We created 0.4—10 keV band light curves by adding the 
MOS1 and MOS2 data only. Although the PN light curves 
have a better signal to noise, they are shorter than the MOS 
light curves. Since the amplitude of the X-ray variations 
in AGN generally increases with increasing time scale, we 
decided to use the longest possible light curvefl We binned 
the combined MOS light curves to intervals of size 500 s. 
In this way the number of photons in each bin is large 
enough (even in the case of the faintest source, NGC 4639) 
to guarantee the applicability of the traditional \ 2 test to 
examine whether a source is variable or not. We accept that 
a light curve shows significant variations if the probability 
of the "null hypothesis" (i.e. the source is not variable) is 
less than 5%. 

The 0.4-10keV, MOS1+2 light curves are shown in 
Fig. EJ NGC 4235 and NGC 4639 do not show signifi- 
cant variations: xl 2 35 — 15.07/21 degrees of freedom (dof), 
and X4 63 g = 26.1/27 dof, respectively. On the other hand 
we do observe significant variations in the NGC 5503 and 
NGC 5273 light curves (xi 033 = 44.3/ 22 and xi 27 3 = 
552/30 dof, respectively). 

We use the normalized excess variance, a^xs ( e -S- 
Nandra et al. 1997) as a measure of the intrinsic variability 



In any case, the PN and MOS light curves look similar in 
the time intervals where both instruments were in operation. 
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Object 


Observing date 


Instrument 


Obsld 


Mode 


Filter 


Exposure 


Count rate 




UT 










(ksec) 


cts/s'-*' 


NGC 4235 


Jun. 09, 2004 


PN 


0204650201 


Large Window 


Thin 


9.6 


0.96 




09:04-12:43 


MOS1 




Prime Partial 


Thin 


10.9 


0.31 






MOS2 




Prime Partial 


Thin 


10.9 


0.31 


NGC 4639 


Dec. 12, 2001 


PN 


0112551001 


Extd. Full Frame 


Medium 


8.8 


0.28 




16:19-18:19 


MOS1 




Full Window 


Thin 


14.2 


0.09 






MOS2 




Full Window 


Thin 


14.2 


0.09 


NGC 5033 


Dec. 18, 2002 


PN 


0094360501 


Full Frame 


Medium 


9.1 


2.40 




15:05-19:01 


MOS1 




Full Window 


Medium 


11.5 


0.68 






MOS2 




Full Window 


Medium 


11.5 


0.68 


NGC 5273 


Jun. 14, 2002 


PN 


0112551701 


Extd. Full Frame 


Medium 


9.3 


1.57 




13:15-16:37 


MOS1 




Full Window 


Thin 


14.2 


0.57 






MOS2 




Full Window 


Thin 


14.2 


0.56 



: Background subtracted values. 



amplitude of the light curves which show significant varia- 
tions. This estimator is defined as, 



'NXS 



S 2 - < aj rr > 
< x > 2 : 



(1) 



where < x >, and S 2 = (1/N data ) £ilT a ( x i - < x >) 2 are 
the mean and variance of the light curve, while < <7 2 rl . >= 

a err i i s ^ ne average contribution of the 
Poisson noise process to the observed scatter around the 
mean (Ndata is the number of points in the light curve) . Its 
square root (the so called "fractional root mean square am- 
plitude" , /rms) indicates the average variability amplitude 
of a source as a fraction of the light curve mean. We find 
that / rms ,5033 = 4 ± 1% and / rms , 5 273 = 18 ± 1% (errors 
account only for the measurement errors in the light curve 
points and have been estimated according to the prescrip- 
tion of Vaughan et al., 2003). 

In order to investigate further the variability of 
NGC 4235 and NGC 4639, we considered ASCA data as 
well. The main advantage of the ASCA light curves is that 
they are are longer than the XMM-Newton light curves. 
In the case of NGC 4235, we employed its 0.5 - 10 keV, 
5760 s binned, SIS 0+1 combined light curves from the 
58 ksec long, December 1998 observation (sequence num- 
ber 76005000; the light curves were downloaded from the 
TARTARUS database). The background subtracted light 
curve is plotted in the upper panel of Fig.O We detect sig- 
nificant variations (x 2 = 39.8/10 dof) with / rms q 5-iokcV = 
6.2 ±1.8%. 

NGC 4639 has also been observed by ASCA. However, 
its background subtracted, 0.5-10 keV, SIS 0+1 combined 
TARTARUS light curve (shown in the bottom panel of 
Fig. [2]) does not reveal any significant variations during its 
~ 40 ksec observation in December 1998 (sequence number 
76007000). 

The top panel in Fig. [3] shows the logarithm of cr^xs 
plotted as a function of the logarithm of the X-ray (i.e., 
2 — 10 keV) luminosity for nearby, bright Seyfert 1 nuclei 
(open circles; data are taken from O'Neill et al. 2005). The 
solid line is the best linear fit to the data. Filled squares 
in the same figure indicate the [log(cr§ rxs ), log(Lx)] data 



0.6 

0.5 



CD 
to 



~ 0.2 



& 0.1 



o 

o 



CD 




1.2 

0.9 
1.5 
1.2 
0.9 
0.6 



, , , | 
NGC 4235 



H 1 — I h 



H 1 1 — h 



H 1 1 — h 



NGC 4639 



-I — I — I — h 



-I — I — I — h 



-I — I — I — h 



-I- 



NGC 5033 



1 — 1 — 1 — 1 — 1 — i 
♦ ♦*/♦ ♦♦♦ 



-I — I — I — h 



H — 1 — 1 — h 
NGC 5273 



-+- 



5000 10000 

Time (sec) 



15000 



Fig. 1. Background subtracted, combined MOS1+2, 0.4- 
10 keV light curves of the four sources, binned in 500 s 
intervals. Time is measured from the start of each obser- 
vation. Filled squares in all panels indicate the respective 
background light curves, rescaled appropriately in each case 
(errors are smaller than the symbols size). 



of NGC 5273, NGC 5033 and NGC 4235 (using the X-ray 
luminosity estimates listed in Table 3). The thick arrow in 
this figure indicates the upper 3cx upper limit of NGC 4639, 
as estimated using the MOS1+2 light curves. The plot sug- 
gests that, apart from NGC 5273, the other three Seyferts 
are less variable than what we would expect, if we consider 
the best line fit to the O'Neill et al. data. 

However, there has been growing evidence the last few 
years that the excess variance in AGN is primarily related 
to BH mass (Papadakis 2004, O'Neill et al. 2005). If that 
is the case, then a "variability amplitude vs luminosity" 
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Fig. 2. Background subtracted, combined SIS 0+1, 0.5 — 10 
keV band light curves of NGC 4235 and NGC 4639, binned 
in 5760 s intervals. Time is measured from the start of each 
observation. 



relation is expected only if the objects in the sample are 
accreting at a similar rate (Papadakis 2004) . In the bottom 
panel of Fig. [31 we plot the logarithm of cr^xs as a function 
of the central black hole mass, Mbh (in Solar mass units), 
for the O'Neill et al. objects and the four Seyferts (BH mass 
estimates are listed in Table 2). This plot shows clearly 
that their variability amplitudes are in agreement with the 
amplitudes we measure in other Seyferts. 

The fact that NGC 5273 shows the largest variability 
amplitude among the four objects in our sample is in agree- 
ment with its small black hole mass. It is not surprising that 
the NGC 4235 MOS light curve does not exhibit significant 
variations, because it is the shortest of the light curves we 
study, and this object has the largest Mbh as well. The 
variations we observe in the ~ 5.5 times longer ASCA light 
curve are of the "correct" amplitude. As for NGC 4639, it is 
rather surprising (given its small black hole mass) that nei- 
ther the MOS nor the ASCA light curves reveal any signifi- 
cant variations. However, this is the faintest among the four 
sources, so the uncertainty in the measurement of its cr NXS 
is quite large. We need a longer XMM-Newtonobseiva.tion 
to detect long term variations in this object. 

Note that, in order to achieve the highest possible 
signal-to-noise ratio, we have used full band light curves 
(i.e. either 0.4 — 10 or 0.5 — 10 keV) to estimate cr^ xs , while 
O'Neill et al. use 2 — 10 keV band light curves. However, this 
should not affect our results significantly, as cf^xs 2-iokeV 
and cr NXS 5 _ 2kcV are comparable in AGN (Nandra et al., 
1997). On the other hand, the NGC 5053 and NGC 5273 
light curves are half the length of the light curves that 
O'Neill et al. use. We therefore conclude that the X-ray 
variability mechanism in these Seyferts operates in the same 
way as in more luminous AGN. 
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Fig. 3. Top panel: Log of cr NXS plotted as a function of log- 
arithm of X-ray luminosity for 33 AGN (open circles; data 
are taken from O'Neill et al. 2005). Filled diamonds indi- 
cate the [log(cr£ xs ), log(Lx)] data of NGC 5273, NGC 5033 
and NGC 4235, while the thick arrow indicates the 3<r up- 
per limit on log(cr Nxg ) for NGC 4639. Bottom panel: The 
log of ct nxs plotted as a function of log(BH mass) (in Solar 
mass units) relation for the 33 AGN of O'Neill et al. and 
the four objects in our sample. 



Table 2. Distance and BH mass estimates. 



Object 


Distance" 


M BH 


ref 6 




(Mpc) 


(M s ) 




NGC 4235 


38.0 


4.8 x 10 7 


(1) 


NGC 4639 


18.4 


7.1 x 10 6 


(2) 


NGC 5033 


15.2 


2.8 x 10 7 


(3) 


NGC 5273 


17.5 


1.6 x 10 6 


(4) 



a Luminosity distances from NED. 

b References for the black hole masses: (1) Dong & de Robertis 
(2006); (2) Panessa et al. (2006); (3) Chiaberge et al. (2005); 
(4) Wu & Han (2001). 



2.3. X-ray Spectral Analysis 

For the X-ray spectral analysis, source spectra were 
grouped with a minimum of 20 counts per energy bin. 
Spectral fits have been performed with the XSPEC vll.3 
package. Spectral responses and the effective area for the 
PN and MOS spectra were generated with the SAS com- 
mands rmfgen and arfgen. Individual fits to the PN and 
the MOS data resulted in similar model parameter values. 
For this reason, we present the results from the joint model 
fits to the PN and the MOS data in the 0.4-10keV energy 
band. In all model fits we add Galactic interstellar absorp- 
tion, using the Nh values reported by Dickey & Lockman 
(1990). The best fitting model parameter values are pre- 
sented together with their 90% errors in Table [3) In the 
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same table we also list our derived 2—10 keV fluxes and 
luminosities. As K-corrections can be neglected for these 
nearby objects we used the luminosity distances as given 
in NED (and listed in Table [2]) for the calculation of the 
luminosities. 



2.3.1. NGC 4639 and NGC 5033 

Cappi et al. (2006) have already studied the XMM-Newton 
observations of these two sources, and our results are in 
agreement with theirs. In the case of NGC 4639, a simple 
power-law (PL) model with T ^ 1.9 plus cold absorption 
(with Nh fixed at the Galactic value) resulted in an accept- 
able fit to the full band spectra of all three instruments 
(Xred/dcL = 1.11/230; best fitting model values are listed 
in Table [3|). We do not find an indication of a cold absorbing 
column in excess of the Galactic value, and the addition of 
a Gaussian line at ~ 6.4 keV does not improve significantly 
the goodness of the fit. 

The full band spectrum of NGC 5033 can be fitted well 
by a power law plus a narrow Gaussian iron line, assuming 
galactic absorption only (the line's width in this 
well as in the other sources below, is kept fixed at 0.01 
keV, appropriate in the case of a narrow line) . The addition 
of the line improves the fit quality by A% 2 = 44.1 for 2 
additional dof. The best model fitting parameters are listed 
in Tabled 



2.3.2. NGC 5273 

The XMM-Newton observation of NGC 5273 has also been 
studied by Cappi et al. (2006). For this object our results 
do not agree with theirs, mainly because they consider the 
spectrum over a smaller energy band. A "PL plus Gaussian 
line plus Galactic absorption" model fits well all the EPIC 
spectra in the 2— lOkeV band. When extrapolated to lower 
energies, a strong, wide absorption feature appears around 
~ 1 keV. A number of complex models suitable to change 
the emission characteristics (like absorption edges, partial 
covering) resulted in fits of moderate quality (x 2 cd — 
which show systematic deviations from the data in certain 
energy ranges. 

Cappi et al. (2006) found that a PL plus a thermal 
plasma model with T ~ 1.4 and kT ~ 0.2 keV can fit well 
the spectra. We find that their model provides a good fit 
only in the 'restricted' energy range of 0.6— lOkeV, while 
at lower energies it diverges dramatically from the data. 
Excellent fits to all spectra over the 0.4— lOkeV energy 
range could, finally, be achieved by a PL plus a Gaussian 
line plus a thermal plasma emission model (parametrized 
with the component MEKAL in XSPEC), modified by the 
presence of an ionized absorber and the cold Galactic mate- 
rial, (i.e. by a WABS xABSORIx (MEKAL + PL + GAUSS) 
model in XSPEC terminology; the ABSORI power-law pho- 
ton index was kept tied to the PL slope). This model yielded 
a x 2 °f 1-05 for 1043 dof. The soft thermal flux in the 
0.5— 2keV band accounts for about 10% of the soft power 
law flux. Fig. [4] shows the above mentioned best model-fit 
to the EPIC data, together with the best fit residuals. 



Table 3. The X-ray model fit results. 



N4235 N4639 N5033 N5273 



N H , Ga i 


1.53 
16 ± 1 


2.35 


1 


0.96 

94 ±4 


PL 

r 

A* 


1.61 
(±0.03) 

5.9 
(±0.2) 


1.86 
(±0.04) 

1.33 
(±0.04) 


1.70 
(±0.01) 

10.3 
(±0.1) 


1.44 
(±0.03) 
10.5+r 5 


GAUSS 

E ; [keV] 

EW[eV] 


6.36 
(±0.02) 

242 +124 


- 


6.43 
(±0.02) 
27511* 


6.40 
(±0.02) 
1901™ 


ABSORI 

e 








26±6 


MEKAL 

kT[keV] 








0.18 
±0.03 


xLi/dof 

rb 

r 2-10keV 
L^-lOkcV 


1.08/596 
2.8 
4.8 


1.11/230 
42.7 
0.17 


1.01/901 
4.3 
1.2 


1.06/1042 
6.3 
2.3 



Model fits were performed to the joint PN ± MOSl±2 spectra 
in the 0.4— lOkeV range. Errors are at the 90% level for one 
interesting parameter. 
f : in units of 10 20 cm" 2 . 

*: PL normalization at 1 keV in 10 -4 photons cm -2 s" 1 keV -1 . 
a The ionization parameter which determines the absorber 
ionization state, and is defined as £ = L/(nR 2 ). 
b in units of 10 _12 ergcm -2 s -1 . 
c in units of 10 41 erg s _1 . 



2.3.3. NGC 4235 

The XMM-Newton observation of this source has not been 
studied in the past. A simple PL model with T ~ 1.6 plus 
a ~ 6.4 Gaussian line fits well the two MOS and the PN 
spectra in the 2 — 10 keV band. Fig. [5] shows the spectra 
together with the 2 — 10 keV best fitting power law model 
extrapolated to low energies, assuming Galactic absorption 
only (upper panel) and the respective best model fitting 
residuals (lower panel) . The flux deficit at energies below ~ 
1 keV clearly indicates the presence of intrinsic absorption 
in this object. 

A PL model plus a Gaussian Line plus a cold absorp- 
tion component (WABS in XSPEC) with Nh left as a free 
parameter gives acceptable fits to the full band, EPIC spec- 
tra of the source. The best model fitting results are listed 
in Table 03 The inclusion of a Gaussian line improves the 
simple PL model fit substantially in the case of the PN 
spectrum only (A% 2 = 22.3 for 2 additional parameters). 
The iron line is not required in the case of the MOS data, 
because of their lower statistical quality. 
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Fig. 4. Simultaneous model fit to the PN, MOS1 (open 
squares) and MOS2 (asterisks) fit to the NGC 5273 0.4- 
10 keV data. The model shown here suggests that the flux 
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Fig. 5. The 0.4-10 keV, EPIC spectra of NGC 4235, plot- 
ted together with the 2 — 10 keV, best model fitting "power- 
law plus Gaussian line" model (upper panel) . The model is 
extrapolated to lower energies assuming Galactic absorp- 
tion only. The respective model fitting residuals plot is 
shown in the lower panel of the figure (symbols are as in 
Fig. H 



3. The near-IR to X-ray SEDs 

In order to construct broad band SEDs we use the results 
from the spectral analysis of the XMM-Newton data, nu- 
clear flux measurements from high-resolution HST images 
(Section 3.1), as well as near-IR measurements from the 
literature. 



Table 4. Near-IR and optical-UV fluxes. 
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2.86xl0~ 16 
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5487A 
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3004A 


6.73xl0- 16 
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1 in erg cm -2 s _1 A~ 1 . 

2 in erg s _1 . 

References: (1) Alonso-Herrero et al. (2003); (2) This work. 

For the near-IR data, we use the nuclear flux mea- 
surements of Alonso-Herrero et al. (2003), derived from 
high-resolution iJSTNICMOS images at 1.1 fim (F110W 
filter) and 1.6 /im (F160W filter). We also use their nu- 
clear flux measurements in the K (A ccn t ra i = 2.12 /im) and 
L-band (A cen t ra i = 3.51 /itm), derived from ~ arcsec resolu- 
tion images taken by the 3.8 m NASA IRTF telescope on 
Mauna Kea. The near-IR data are summarized in Table [U 
Luminosities at each waveband have been estimated as de- 
scribed in Section 2.3. 

Finally, we searched the MAST archive at STScI for all 
available WFPC2 images for the four galaxies, in order to 
measure their nuclear fluxes in as many as possible optical 
wavebands. We describe below the HST data selection and 
reduction methods we have used. 

3.1. The HST Observations 

All four Seyferts have been observed in at least one fil- 
ter with WFPC2 on board HST. A summary of the HST 
observations we used is given in Table O The date of obser- 
vation, exposure time and HST proposal number are listed 
in Columns 2, 4 and 5. Column 3 lists the WFPC2 filters, 
together with the respective A cen trai (in parenthesis). For 
the galaxies with more than one image in each filter, we 
chose to study the one with the largest exposure which did 
not cause saturation effects in the central region. 

The WFPC2 images were already processed by the stan- 
dard STScI reduction pipeline (Baggett et al. 2002). Since 
we have chosen images with no saturation effects, the only 



Papadakis et al.: Low lum 
Table 5. Summary of the HST Observations. 



Object 


Observation 
Date 


Filter 


Exp. 
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Prop. 
ID # 


NGC4235 


04/29/1995 


F606W(6030A) 


500 


5479 


NGC4639 


01/15/1995 
01/15/1995 
01/15/1995 


F547M(5487A) 
F439W(4316A) 
F300W(3004A) 


230 
700 
1000 


5381 
5381 
5381 


NGC5033 


05/28/2001 
12/18/1994 
12/18/1994 


F606W(6030A) 
b 547M(5487A) 
F439W(4316A) 


400 
230 
700 


8597 
5381 
5381 


NGC5273 


04/06/1997 
07/04/2000 
01/14/1995 
01/14/1995 
01/13/1995 


F791W(7898A) 
F606W(6030A) 
F547M(5487A) 
F439W(4316A) 
F300W(3004A) 


40 
400 
230 
700 
1000 


6419 
8597 
5381 
5381 
5381 



additional processing step required is the removal of cos- 
mic rays. For that reason, we used the filter/cosmic task of 
the MIDAS astronomical package. Any residual cosmic ray 
events as well as bright foreground stars were removed by 
hand. 

In order to measure nuclear fluxes we derived the radial 
brightness profile of each galaxy. To this end, we performed 
ellipse fitting to the isophotes of the galaxy images. This 
choice is motivated by the fact that the isophotes of galax- 
ies, especially elliptical (E) and lenticular (SO) as well as 
the bulge of spiral galaxies, are not far from ellipses. This 
technique has been widely used in the past by various au- 
thors, mainly as a method of retrieving embedded galaxy 
structures that are hidden by the large-scale distribution 
of light of the main body of the galaxy (see Xilouris & 
Papadakis 2004, and references therein). Using the fit/ell3 
task of MIDAS, which is based on the formulas of Bender 
& Mollenhoff (1987), we fitted the isophotes of the galaxy 
with ellipses and thus derived the average surface brightness 
profile along the major axis of the galaxy. To convert from 
counts to physical units, we multiplied the pixel counts by 
the value of the keyword PHOTFLAM in the image headers 
and divide by the value of the EXPTIME keyword. 

We also created synthetic PSF's using the Tiny Tim 
software (Krist 1995), which, for WFPC2, produces an 
accurate representation of the central region of the PSF, 
thus appropriate for our purpose. We aligned the synthetic 
PSF's with the position of the nucleus, and constructed the 
radial profile of the PSF's with the same method as we do 
for the galaxies. 

Since we are not interested in modeling the bulge of the 
galaxies on large scales, we considered the radial bright- 
ness profile only in the central 2" and we fitted it with a 
"Sersic's function + PSF" model, i.e. a model of the form 
I(R) = 7(0) exp{-(i?/i? ) 1/n } + K x PSF(i?), where I(R) 
is the surface brightness, 1(0), Rq, and n are constants to be 
determined during the fitting process, PSF(i?) is the PSF 
radial profile, and if is a constant which in effect determines 
the unresolved nuclear flux. Fig. [5] shows as an example the 
brightness profile of the innermost 5" in NGC 4235, to- 
gether with the best model fitting line. 



le-14 



le-15 



le-16 



NGC 4235 




Fig. 6. Radial brightness profile of the central arcsec of 
NGC 4235 (from the F6006W image). The solid line shows 
the best fit to the profile (the fit was performed to the 
innermost 2" only)). 



Our results are listed in Table 2] The nuclear fluxes 
are corrected for Galactic extinction, adopting the values 
of Schlegel et al. (1998). In the case of NGC 4235 we also 
considered intrinsic absorption, given the clear presence of 
a dust lane in its HST image and the X-ray spectral fitting 
results which indicate absorption of the nuclear source by 
material within the galaxy. For the reddening correction in 
the optical band we converted the best fitting Nh value to 
A v using the relation A v = N H /2.2 x 10 21 (Ryter 1996), 
and we found the extinction at 6030A using the relations 
of Cardelli et al. (1989). 

Ho & Peng (2001) have also used the WFPC2/F547M 
HST images of NGC 4639, NGC 5033, and NGC 5273 to 
determine their nuclear flux, using an advanced method for 
galaxy image decomposition. They report a nuclear flux 
of 3.9 x 10"17, 9 x 10~16, and 2.7 x 10"16 erg cm" 2 s" 1 , 
respectively. These estimates are in agreement with our re- 
sults. 



3.2. The overall SEDs 

Fig. [7] shows the near IR-to- X-ray SEDs of the four 
Seyferts. The solid line is the average SED of radio-quiet 
quasars from Elvis et al. (1994); dotted and dashed lines 
show the 68 and 90 percentile contours on each side of 
the mean SED, respectively. The average quasar SED is 
normalized to the flux of NGC 4235 in the F606W filter, 
and to the flux of the other three objects in the F439W 
filter. Note that all the measurements plotted in Fig. [7] 
are non-simultaneous. They were taken several years apart, 
and variations (by factors of a few, specially in the X-ray 
band) over these time scales are common in AGN. This fact 
could introduce uncertainties which should be kept in mind 
when comparing the non-simultaneous measurements with 
the mean quasar SED. 

The near IR/optical SEDs of NGC 5033 and NGC 5273 
show a concave spectrum, which is almost identical, in 
shape, to the mean quasar SED. This spectral shape is often 
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Fig. 7. The near-IR through the X-ray SED for the four objects 
in our sample. Upper limits are 3a. The solid, dotted and dashed 
lines indicate the mean radio-quiet quasar SED from Elvis et al., 
and the the 68 and 90 percentile contours on each side. Filled 
squares are our measurements from the HST data, open squares 
are IR data from the literature. Crosses indicate the nuclear 
SED of NGC 3998 (see text fro details). 



interpreted as the signature of the presence of the "stan- 
dard", geometrically thin/optically thick disc emission and 
of dust heated by the central engine. In NGC 4639 we ob- 
serve the optical part of the spectrum only, which is rising 
towards the UV (as expected if we observe thermal emission 
from the accretion disc) . 

In the case of NGC 4235 and NGC 5033, the X-ray 
to optical luminosity ratio is almost identical to what we 
observe in quasars. In NGC 4639 and NGC 5273, the 
optical-UV bump appears to be weaker than the X-rays 
by a factor of a few. Since it is in the X-ray band where 
AGN show their largest amplitude variations, this differ- 
ence may not be intrinsic. For example, during its XMM- 
TVewfonobservation, NGC 4639 appears to be ~ 3 and ~ 10 
times brighter than during its ASCA(Ter&shima, et al. 2002) 
and CHANDRA observation (Ho et al. 1999). 

On the other hand, this difference may be intrinsic, in- 
dicating that accretion in these objects operates in a mode 
which is different than that in powerful quasars. To illus- 
trate this point, the crosses in all panels of Fig. [7J show the 
nuclear SED of NGC 3998. Data were taken from Table 4 
of Ptak et al. (2004) . In the case there were multiple mea- 
surements at a certain wavelength, we chose the ones de- 
rived from the smallest aperture size. The NGC 3998 SED 
is rescaled to the optical luminosity of the objects in our 
sample. Figure [7] shows that the optical to X-ray SEDs of 



NGC 3998, NGC 4235 and NGC 4639 (at least) are almost 
identical. 

The NGC 3998 SED is consistent with RIAF models. 
Ptak et al. (2004) found that, although a "pure" thin-disc 
model does not account the mid-IR emission in NGC 3998, 
a truncated thin-disc (at around ~ 300 Schwarzschild 
radii), surrounding an inner RIAF, fitted its mid-IR to X— 
ray spectrum well (when the accretion rate is ~ 10~ 3 of the 
Eddington limit). Consequently, the qualitative comparison 
of the (rather sparse) SEDs of the objects in our sample 
with either the average SED of the radio-quiet quasars or 
SEDs which are consistent with RIAF models, does not al- 
low us to reach solid conclusions regarding the accretion 
mode in these objects. 

We reach more solid results when we compare the op- 
tical and X-ray luminosity, in quasars and LLAGNs, in a 
more quantitative way. Fig. [5] shows the monochromatic lu- 
minosity Ly at 5400A, plotted as a function of Lx, for PG 
quasars (data taken from Elvis et al., 1994). Filled squares 
show the monochromatic luminosity at 5487A plotted as 
a function of Lx in the case of the fours Seyferts of our 
sample (using the data listed in Tables 3 and 4), while the 
open squares show data for other nearby Seyferts, using the 
nuclear flux 5487A and X-ray flux measurements of Ho & 
Peng (2001) and Cappi et al. (2006). The cross in the same 
figure indicates the 5000A and average 2-10 keV luminos- 
ity of NGC 3998, estimated using the flux measurements of 
Ptak et al. 2004 (as listed in their Table 2) and assuming a 
distance of 14.1 Mpc. 

We find that a line of the form Ly oc l.l(±0.1)Lx (solid 
line in Fig. |5J) fits the quasar data well. The dashed line in 
the same figure shows its extrapolation to lower luminosi- 
ties. The agreement between the low luminosity Seyferts 
and the powerful PG quasars is excellent. We do not ob- 
serve any sign of a "dichotomy" in this plot, over seven 
orders of magnitude in X-ray luminosity, which might indi- 
cate that the low luminosity AGN accrete in the inefficient 
radiative regime. In fact, the addition of the low luminosity 
objects seems to strengthen the "optical vs X-ray luminos- 
ity" correlation we detect among the PG quasars alone. We 
conclude that the optical to Xray SEDs of the four objects 
in our sample (and other nearby, low luminosity Seyferts) 
is not different from the mean SED of quasars, which are 
two to five orders of magnitude more luminous. 

4. Discussion and Conclusions 

We have used archival XMM-Newton&nd ASCAd&ta, in or- 
der to measure the X-ray variability amplitude of four low 
luminosity Seyfert 1-1.5 galaxies, and compare it with the 
amplitude of brighter AGN. We have used the same XMM- 
Newtonda,ta, to determine accurately their nuclear X-ray 
luminosity and spectral shape, and HST data to measure 
the optical/UV flux of their active nucleus. We then com- 
bined the XMM-Newtonspectra., the HST measurements, 
and near-IR data from literature to construct their broad 
band, near-IR to X-ray SED and compare it with the mean 
SED of radio-quiet quasars. Our main aim is to investigate 
whether accretion in these objects operates in a mode which 
is different from the accretion mode in their more luminous 
counterparts. 

Such comparisons are frequently used to investigate the 
nature of the nuclear source in low luminosity AGNs. Most 
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Fig. 8. The Ly versus Lx plot for PG quasars (open cir- 
cles), nearby Seyferts (open squares), the four Seyferts in 
our sample (filled squares), and NGC 3998 (cross). The 
solid line indicates the best linear fit to the quasar sample, 
and the dashed line is its extrapolation to lower luminosi- 
ties. 



of the similar work in the past has focused mainly on the 
comparison between the properties of galaxies with low- 
ionization nuclear emission-line regions (LINERs; Heckman 
1980) and luminous AGN. However, the nature of LINERS 
and their relation to AGN has been debated for several 
years. For this reason, we decided to study objects with 
a firm optical classification as a Seyfert 1-1.5 galaxy. In 
this way, we are certain that they host an accretion driven, 
AGN-like, UV to X-ray nuclear source. In addition, the 
clear presence of broad lines in their optical spectra implies 
that we can view their central source directly, and hence 
determine their nuclear SED reliably. 

Furthermore, using the X-ray luminosities listed in 
Table 3, and the BH mass estimates listed in Table 2, we 
find that the (L x /L E dd) ratio of NGC 5273 is 6 x 1CT 4 . 
The value of the same ratio in the case of NGC 4235, 
NGC 4639 and NGC 5033 is even smaller: 7.7,1.8 and 
3.3 x 10 -5 , respectively. This is much smaller than the 
same ratio of powerful quasars, and even 20 to 80 times 
smaller than the ratio of other Seyfert 1 galaxies in the 
Palomar survey (like for example NGC 4051, NGC 4151 
and NGC 5548). It is also below the limit of 10~ 3 where 
Merloni et al. (2003) suggest that the accretion mode may 
change from a thin, optically thick disc to that of a RIAF. 
We also converted the measured F547 fluxes of NGC 4639, 
NGC 5033 and NGC 5273, and the measured F606W flux 
of NGC 4235, to flux at 7000 A, assuming an optical spec- 
tral index of 1, and we estimated the monochromatic lumi- 
nosity at this wavelength. Chiaberge et al. (2005) suggest 
that if Lo/LEdd < 10~ 4 , then a RIAF may operate in an 
AGN. In the case of NGC 4235, NGC 4639 and NGC 5033, 
this ratio is 6 - 60 x 10 -6 , while in the case of NGC 5273, 
Lo/LEdd ~ 1-5 x 10- 4 . 



In summary, the four objects of our study i) should host 
a relatively unobscured AGN-like nucleus, ii) are of low 
luminosity, and iii) accrete on a low rate as well. Hence 
they are good candidates to detect radiation from RIAF- 
like processes. 

We detect significant X-ray variations in the XMM- 
NewtonligU curves of NGC 5033 and NGC 5273, and in 
the AS C Alight curve of NGC 4235. The variability ampli- 
tude is 5-15% of the light curve's mean, and is exactly of the 
"right" order, given their BH mass. In agreement with Ptak 
et al. (1998), we find that these objects do not follow the 
"variability amplitude - luminosity" relation of other more 
luminous AGN. However, this result does not necessarily 
imply a larger characteristic size for the X-ray producing 
region than is the case in "normal" AGN. On the contrary, 
if the size is determined by BH mass, as expected, then our 
results suggest that, when scaled to BH mass, the size of 
their X-ray sources is similar to the size of the nuclei in 
more luminous AGN. Their "misplacement" in the "ampli- 
tude - luminosity" plane is not due to their variability am- 
plitude being too low but due to their luminosity being low, 
most probably because of their low accretion rate. Their X- 
ray spectral shape is also consistent with the shape of more 
luminous AGN. Even the best fitting spectral slope in the 
case of NGC 5273 is not exceptionally unusual; NGC 4151, 
a well studied, luminous Seyfert galaxy, also shows a very 
flat X-ray spectrum of L ~ 1.5 (see e.g. de Rosa et al. 
2007). We therefore conclude that, from the X-ray point 
of view, the central engine in these low luminosity Seyferts 
operates in the same way as in the more luminous AGN. 

The near-IR to X-ray SEDs of the objects in our sample 
are rather sparse, and their qualitative comparison with ei- 
ther the average SED of radio-quiet quasars or the SED 
of other low-luminosity objects which appear to host a 
RIAF, is not conclusive. If we take into account the fact 
that we used measurements taken several years apart, and 
the possibility of large amplitude variations, at least in 
the X-ray band, then it is reasonable to conclude that the 
SEDs of the AGN in our sample are similar to the average 
SED of quasars. On the other hand, the SED of at least 
NGC 4235 and NGC 4639 are almost identical to the SED 
of NGC 3998, which is consistent with a model of a trun- 
cated disc and a RIAF at smaller radii (see also the case 
of NGC 4579 and M81 whose SEDs are consistent with a 
similar model; Quataert et al. 1999). 

The situation becomes clearer when we compare the op- 
tical luminosity, Ly, with the X-ray luminosity, Lx, in the 
4 Seyferts of our sample (and a few other LLAGNs) and in 
quasars (i.e. AGN which are 10 2 to 10 5 times more lumi- 
nous). We find that the quasar optical luminosity increases 
almost proportionally with Lx- The extrapolation of the 
quasar "Ly vs Lx" best fit line to lower luminosities agrees 
very well with the LLAGN data. It is then reasonable to 
assume that, the observed "continuity" in the (Ly, Lx) dis- 
tribution of quasars and LLAGNs, also suggests that the 
optical/X-ray SEDs in AGN are similar, hence, there is 
no accretion mode "phase transition" , over almost seven 
orders of magnitude. 

Our results arc in agreement with the results of 
Panessa et al. (2006) who study the X-ray and opti- 
cal emission line correlation in a sample of AGN with a 
wide range of Eddington ratios, and conclude that low- 
luminosity Seyferts are simply scaled-down versions of lumi- 
nous quasars. Similar results were obtained by Maoz (2007) 
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who studied the SEDs of 13 low luminosity LINERs, and 
found that radiative efficient discs may operate even in 
these objects, which are characterized by very low accre- 
tion rates. 

Although we have tried to use the best available data to 
determine the SED for the four Seyferts, the fact that the 
measurements we use in the various bands are separated 
by years, introduces some uncertainty. Significant improve- 
ment can be achieved if we manage to determine the average 
nuclear flux of each object in various bands. This can be 
achieved only if we observe individual objects many times 
over a time scale of a few years at least, which is not an 
easy task. However, the "variability amplitude test" results 
are not affected by these limitations. Even if we consider 
only them, the main result of our work is that accretion op- 
erates in the "normal" way even in objects where the ratio 
Lx/LEdd is as low as a fewxlO -5 . 

Perhaps, it is not just the accretion rate that determines 
the accretion mode in AGN. Black hole spin and/or mass 
may also play a role as to how accretion operates in active 
galaxies. We plan to investigate this issue by studying the 
variability properties of a large sample of low luminosity 
AGNs, with a wide range of BH masses. A dichotomy in 
the "variability amplitude - BH mass" plane, could provide 
us with important clues as which are the main parameters 
which dictate the "switch" in the AGN accretion mode. 
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